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Abstract 

It is important to have fast and effective methods for simplifying 3-manifold triangula- 
tions without losing any topological information. In theory this is difficult: we might need 
to make a triangulation super-exponentially more complex before we can make it smaller 
than its original size. Here we present experimental work suggesting that for 3-sphere tri- 
angulations the reality is far different: we never need to add more than two tetrahedra, 
and we never need more than a handful of local modifications. If true in general, these 
extremely surprising results would have significant implications for decision algorithms and 
the study of triangulations in 3-manifold topology. 

The algorithms behind these experiments are interesting in their own right. Key tech- 
niques include the isomorph-free generation of all 3-manifold triangulations of a given size, 
polynomial-time computable signatures that identify triangulations uniquely up to isomor- 
phism, and parallel algorithms for studying finite level sets in the infinite Pachner graph. 

ACM classes F.2.2; G.2.1; G.2.2; D.1.3 

Keywords Triangulations, 3-manifolds, Pachner moves, 3-sphere recognition, isomorph- 
free enumeration 

1 Introduction 

Triangulations of 3-manifolds are ubiquitous in computational knot theory and low-dimensional 
topology. They are easily obtained and offer a natural setting for many important algorithms. 

Computational topologists typically allow triangulations in which the constituent tetrahedra 
may be "bent" or "twisted", and where distinct edges or vertices of the same tetrahedron 
may even be joined together. Such triangulations (sometimes called semi-simplicial or pseudo- 
triangulations) can describe rich topological structures using remarkably few tetrahedra. For 
example, the 3-dimensional sphere can be built from just one tetrahedron, and more complex 
spaces such as non-trivial surface bundles can be built from as few as six 15_. 

An important class of triangulations is the one-vertex triangulations, in which all vertices 
of all tetrahedra are joined together as a single point. These are simple to obtain [TTJIIl], and 
they are often easier to deal with both theoretically and computationally [SJ [TUl [H] . 

Keeping the number of tetrahedra small is crucial in computational topology, since many 
important algorithms are exponential (or even super-exponential) in the number of tetrahedra 
[51 [^. To this end, topologists have developed a rich suite of local simplification moves that 
allow us to reduce the number of tetrahedra without losing any topological information [31 |TB] . 
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The most basic of these are the four Pachner moves (also known as bistellar moves). These 
include the 3-2 move (which reduces the number of tetrahedra but preserves the number of 
vertices), the 4-1 move (which reduces both numbers), and also their inverses, the 2-3 and 1-4 
moves. It is known that any two triangulations of the same closed 3-manifold are related by a 
sequence of Pachner moves |21) . Moreover, if both are one- vertex triangulations then we can 
relate them using 2-3 and 3-2 moves alone [H]. 

However, little is known about how difficult it is to relate two triangulations by a sequence 
of Pachner moves. In a series of papers, Mijatovic develops upper bounds on the number 
of moves required for various classes of 3-manifolds [171 Ull [HI ISO] . All of these bounds are 
super-exponential in the number of tetrahedra, and some even involve exponential towers of 
exponential functions. For relating one-vertex triangulations using only 2-3 and 3-2 moves, no 
explicit bounds are known at all. 

In this paper we focus on one-vertex triangulations of the 3-sphere. Here simplification is 
tightly linked to the important problem of 3-sphere recognition, where we are given an input 
triangulation T and asked whether T represents the 3-sphere. This problem plays an key role 
in other important topological algorithms, such as connected sum decomposition |11[ 112) and 
unknot recognition [5], and it is now becoming important in computational 4^-manifold topology. 
We can use Pachner moves for 3-sphere recognition in two ways: 

• They give us a direct 3-sphere recognition algorithm: try all possible sequences of Pachner 
moves on T up to Mijatovic's upper bound, and return "true" if and only if we reach one 
of the well-known "canonical" 3-sphcrc triangulations with one or two tetrahedra. 

• They also allow a hybrid recognition algorithm: begin with a fast and/or greedy procedure 
to simplify T as far as possible within a limited number of moves. If we reach a canon- 
ical 3-sphere triangulation then return "true" ; otherwise run a more traditional 3-sphere 
recognition algorithm on our new (and hopefully simpler) triangulation. 

The direct algorithm lies well outside the realm of feasibility: Mijatovic's bound is super- 
exponential in the number of tetrahedra, and the running time is at least exponential in Mija- 
tovic's bound. Current implementations [3] use the hybrid method, which is extremely effective 
in practice. Experience suggests that when T is the 3-sphere, the greedy simplification almost 
always gives a canonical triangulation. If simplification fails, we revert to the traditional algo- 
rithm of Rubinstein |22| : although this runs in exponential time, recent improvements by several 
authors have made it extremely effective for moderate-sized problems [i [a [III [23] 

Our aims in this paper are: 

• to measure how easy or difficult it is in practice to relate two triangulations of the 3-sphere 
using Pachner moves, or to simplify a 3-sphere triangulation to use fewer tetrahedra; 

• to understand why greedy simplification techniques work so well in practice, despite the 
prohibitive theoretical bounds of Mijatovic; 

• to investigate the possibility that Pachner moves could be used as the basis for a direct 
3-sphere recognition algorithm that runs in sub-exponential or even polynomial time. 

Fundamentally this is an experimental paper (though the theoretical underpinnings are in- 
teresting in their own right). Based on an exhaustive study of ~ 150 million triangulations 
(including ^ 31 million one- vertex triangulations of the 3-sphere), the answers appear to be: 

^See [7] for explicit measurements of running time. 
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• we can relate and simplify one- vertex triangulations of the 3-spliere using remarkably few 
Pachner moves; 

• both procedures require us to add at most two extra tetrahedra, which explains why greedy 
simplification works so well; 

• the number of moves required to simplify such a triangulation could also be bounded by 
a constant, which means polynomial-time 3-sphere recognition may indeed be possible. 

These observations are extremely surprising, especially in light of Mijatovic's bounds. If 
they can be proven in general — yielding a polynomial-time 3-sphere recognition algorithm — this 
would be a significant breakthrough in computational topology. 

In Section [5] we outline preliminary concepts and introduce the Pachner graph, an infinite 
graph whose nodes represent triangulations and whose arcs represent Pachner moves. This 
graph is the framework on which we build the rest of the paper. We define simplification paths 
through the graph, as well as the key quantities of length and excess height that we seek to 
measure. 

We follow in Section |3] with two key tools for studying the Pachner graph: an isomorph-free 
census of all closed 3-manifold triangulations with < 9 tetrahedra (which gives us the nodes of 
the graph), and isomorphism signatures of triangulations that can be computed in polynomial 
time (which allow us to construct the arcs of the graph). 

Section 0] describes parallel algorithms for bounding both the length and excess height of 
simplification paths, and presents the highly unexpected experimental results outlined above. 
We finish in Section [5] with a discussion of the implications and consequences of these results. 

2 Triangulations and the Pachner graph 

A 3-manifold triangulation of size rt is a collection of n tetrahedra whose 4n faces are affinely 
identified (or "glued together") in 2n pairs so that the resulting topological space is a closed 
3-manifoldll We are not interested in the shapes or sizes of tetrahedra (since these do not affect 
the topology), but merely the combinatorics of how the faces are glued together. Throughout 
this paper, all triangulations and 3-manifolds are assumed to be connected. 

We do allow two faces of the same tetrahedron to be identified, and we also note that distinct 
edges or vertices of the same tetrahedron might become identified as a by-product of the face 
gluings. A set of tetrahedron vertices that are identified together is collectively referred to as a 
vertex of the triangulation] we define an edge or face of the triangulation in a similar fashion. 

Figure [1] illustrates a 3-manifold triangulation of size n — 2. Here the back two faces of 
the first tetrahedron are identified with a twist, the front faces of the first tetrahedron are 
identified with the front faces of the second using more twists, and the back faces of the second 
tetrahedron are identified together by directly "folding" one onto the other. This is a one-vertex 
triangulation since all eight tetrahedron vertices become identified together. The triangulation 
has three distinct edges, indicated in the diagram by three distinct arrowheads. 

Mijatovic [17j describes a canonical triangulation of the 3-sphere of size n — 2, formed by a 
direct identification of the boundaries of two tetrahedra. In other words, given two tetrahedra 

■^It is sometimes useful to consider bounded triangulations where some faces are left unidentified, or ideal 
triangulations where the overall space only becomes a 3-manifold when we delete the vertices of each tetrahedron. 
Such triangulations do not concern us here. 
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Figure 1 : A 3-nianifold triangulation of size n — 2 



ABCD and A'B'C'D', we directly identify face ABC with A'B'C, ABD with A'B'D', and so 
on. The resulting triangulation has four faces, six edges, and four vertices. 

The four Pachner moves describe local modifications to a triangulation. These include: 

• the 2-3 move, where we replace two distinct tetrahedra joined along a common face with 
three distinct tetrahedra joined along a common edge; 

• the 1-4 move, where we replace a single tetrahedron with four distinct tetrahedra meeting 
at a common internal vertex; 

• the 3-2 and 4-1 moves, which are inverse to the 2-3 and 1-4 moves. 



These four moves are illustrated in Figure [2] Essentially, the 1-4 and 4-1 moves retriangulate 
the interior of a pyramid, and the 2-3 and 3-2 moves retriangulate the interior of a bipyramid. It 
is clear that Pachner moves do not change the topology of the triangulation (i.e., the underlying 
3-manifold remains the same). Another important observation is that the 2-3 and 3-2 moves do 
not change the number of vertices in the triangulation. 

Two triangulations are isomorphic if they are identical up to a relabelling of tetrahedra and a 
reordering of the four vertices of each tetrahedron (that is, isomorphic in the usual combinatorial 
sense). Up to isomorphism, there are finitely many distinct triangulations of any given size. 

Pachner originally showed that any two triangulations of the same closed 3-manifold can be 
made isomorphic by performing a sequence of Pachner moves [21| lf| Matveev later strengthened 
this result to show that any two one-vertex triangulations of the same closed 3-manifold with at 
least two tetrahedra can be made isomorphic through a sequence of 2-3 and/or 3-2 moves [14]. 

^As Mijatovic notes, Pachner's original result was proven only for true simplicial complexes, but it is easily 
extended to the more flexible definition of a triangulation that we use here 1171 . The key step is to remove 
irregularities by performing a second barycentric subdivision using Pachner moves. 




(a) The 2-3 and 3-2 moves 



(b) The 1-4 and 4-1 moves 



Figure 2: The four Pachner moves for a 3-manifold triangulation 
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The two-tetrahedron condition is required because it is impossible to perform a 2-3 or 3-2 move 
upon a one-tetrahedron triangulation (each move requires two or three distinct tetrahedra) . 

In this paper we introduce the Pachner graph, which describes how distinct triangulations of 
a closed 3-manifold can be related via Pachner moves. We define this graph in terms of nodes and 
arcs, to avoid confusion with the vertices and edges that appear in 3-manifold triangulations. 

Definition (Pachner graph). Let M be any closed 3-manifold. The Pachner graph of M, 
denoted ^(M), is an infinite graph constructed as follows. The nodes of ^(M) correspond to 
isomorphism classes of triangulations of M. Two nodes of ^(M) are joined by an arc if and 
only if there is some Pachner move that converts one class of triangulations into the other. 

The restricted Pachner graph of M, denoted 3^i{M), is the subgraph of ^(M) defined by 
only those nodes corresponding to one-vertex triangulations. The nodes of 3^{M) and ^i{M) 
are partitioned into finite levels 1, 2,3,.. ., where each level n contains the nodes corresponding 
to n-tetrahedron triangulations. 
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Level 3 
Level 2 

• Level 1 

Figure 3: Levels 1-3 of the restricted Pachner graph of the 3-sphere 

It is clear that the arcs are well-defined (since Pachner moves are preserved under isomor- 
phism), and that arcs do not need to be directed (since each 2-3 or 1-4 move has a corresponding 
inverse 3-2 or 4-1 move). In the full Pachner graph ^(M), each arc runs from some level i to 
a nearby level i ± 1 or i ± 3. In the restricted Pachner graph ^i(M), each arc must describe a 
2-3 or 3-2 move, and must run from some level i to an adjacent level i ± 1. Figure [3] shows the 
first few levels of the restricted Pachner graph of the 3-sphere. 

We can now reformulate the results of Pachner and Matveev as follows: 

Theorem 1 (Pachner, Matveev). The Pachner graph of any closed 3-manifold is connected. If 
we delete level 1, the restricted Pachner graph of any closed 3-manifold is also connected. 

To simplify a triangulation we essentially follow a path through 3^{M) or ^i(M) from a 
higher level to a lower level, which motivates the following definition. 

Definition (Simplification path). A simplification path is a directed path through either 3^{M) 
or ^i(Af) from a node at some level « to a node at some lower level < i. The length of a 
simplification path is the number of arcs it contains. The excess height of a simplification path 
is the smallest /i > for which the entire path stays in or below level i + h. 

Both the length and excess height measure how difficult it is to simplify a triangulation: the 
length measures the number of Pachner moves, and the excess height measures the number of 
extra tetrahedra required. For the 3-sphere, the only known bounds on these quantities are the 
following: 
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Theorem 2 (Mijatovic |17p. Any triangulation of the 3-sphere can be converted into the canon- 
ical triangulation using less than 6 • 10^n^2^'^'' " Pachner moves. 

Corollary 3. In the Pachner graph of the 3-sphere, from any node at level n > 2 there is a 
simplification path of length less than 6-10^n^2^'^*^ " and excess height less than 3-10^n^2^'^'^ " . 

In the restricted Pachner graph, no exphcit bounds on these quantities are known at all. 



3 Key tools 

Experimental studies of the Pachner graph are difficult: the graph itself is infinite, and even 
the finite level sets grow super-exponentially in size. By working with isomorphism classes of 
triangulations, we keep the level sets considerably smaller than if we had used labelled trian- 
gulations instead. However, the trade-off is that both the nodes and the arcs of the graph are 
more difficult to construct. 

In this section we outline two key algorithmic tools for studying the Pachner graph: a census 
of triangulations (which enumerates the nodes at each level), and polynomial-time computable 
isomorphism signatures (which allow us to construct the arcs). 



3.1 A census of triangulations 

To enumerate the nodes of Pachner graphs, we build a census of all 3-manifold triangulations of 
size n < 9, with each triangulation included precisely once up to isomorphism. Because we are 
particularly interested in one- vertex triangulations as well as triangulations of the 3-sphere, we 
extract such triangulations into separate censuses with the help of the highly optimised 3-sphere 
recognition algorithm described in [7] . The final counts are summarised in Table [1] 



Number of 


All closed 3-manifolds 


3-spheres only 


tetrahedra 


No constraints 


1-vertex only 


No constraints 


1-vertex only 


1 


4 


3 


2 


1 


2 


17 


12 


6 


3 


3 


81 


63 


32 


20 


4 


577 


433 


198 


128 


5 


5184 


3 961 


1903 


1297 


6 


57 753 


43 584 


19 935 


13 660 


7 


722 765 


538 409 


247 644 


169 077 


8 


9 787 509 


7 148 483 


3 185 275 


2142 197 


9 


139 103 032 


99 450 500 


43 461431 


28 691 150 


Total 


149 676 922 


107 185 448 


46 916 426 


31017 533 



Table 1: Counts for 3-manifold triangulations of various types in the census 



The algorithms behind this census are sophisticated; see [3] for some of the techniques 
involved. The constraint that the triangulation must represent a 3-manifold is critical: if we 
just enumerate all pairwise identifications of faces up to isomorphism, there are at least 

[(4n - 1) X (4n - 3) X • • • X 3 X 1] • 6^" o i le 



24" 



2.35 X 10' 
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possibilities for n 9. To enforce the 3- manifold constraint we use a modified union-find 
algorithm that tracks partially-constructed edge links and vertex links; see [3] for details. 
Even with this constraint, we can prove that the census grows at a super-exponential rate: 

Theorem 4. The number of distinct isomorphism classes of 3-manifold triangulations of size 
n grows at an asymptotic rate o/exp(0(nlogn)). 

The proof is detailed, and is given in the appendix. 

For the largest case n = 9, the enumeration of all 3-manifold triangulations up to isomor- 
phism required ^ 85 days of CPU time as measured on a single 1.7 GHz IBM PowcrS processor 
(though in reality this was reduced to 2-3 days of wall time using 32 CPUs in parallel) . The time 
required to extract all 3-sphere triangulations from this census was negligible in comparison. 

3.2 Isomorphism signatures 

To construct arcs of the Pachner graph, we begin at a node — that is, a 3-manifold triangulation 
T — and perform Pachner moves. Each Pachner move results in a new triangulation T', and our 
main difficulty is in deciding which node of the Pachner graph represents T'. 

A naive approach might be to search through nodes at the appropriate level of the Pachner 
graph and test each corresponding triangulation for isomorphism with T'. However, this ap- 
proach is infeasible: although isomorphism testing is fast (as we prove below), the sheer number 
of nodes at level n of the graph is too large (see Theorem HJ . 

What we need is a property of the triangulation T' that is easy to compute, and that uniquely 
defines the isomorphism class of T'. This property could be used as the key in a data structure 
with fast insertion and fast lookup (such as a hash table or a red-black tree), and by computing 
this property we could quickly jump to the relevant node of the Pachner graph. 

Here we define such a property, which we call the isomorphism signature of a triangulation. 
In Theorem [S] we show that isomorphism signatures do indeed uniquely define isomorphism 
classes, and in Theorem [7] we show that they are small to store and fast to compute. 

A labelling of a triangulation of size n involves: (i) numbering its tetrahedra from 1 to n 
inclusive, and (ii) numbering the four vertices of each tetrahedron from 1 to 4 inclusive. We also 
label the four faces of each tetrahedron from 1 to 4 inclusive so that face i is opposite vertex i. 
A key ingredient of isomorphism signatures is canonical labellings, which we define as follows. 

Definition (Canonical labelling). Given a labelling of a triangulation of size n, let Atj denote 
the tetrahedron which is glued to face / of tetrahedron t (so that Atj £ {!,..., n} for all 
t = 1, . . . , n and / = 1, . . . , 4). The labelling is canonical if, when we write out the sequence 
^1,2, ^i,3i ^i,4j ^2,1, • • ■ : ^n,4j the foUowiug properties hold: 

(i) For each 2 < i < j, tetrahedron i first appears before tetrahedron j first appears. 

(ii) For each i > 2, suppose tetrahedron i first appears as the entry Atj. Then the correspond- 
ing gluing uses the identity map: face / of tetrahedron t is glued to face / of tetrahedron 
i so that vertex v of tetrahedron t maps to vertex v of tetrahedron i for each v ^ f. 

As an example, consider the triangulation of size n = 3 described by Table H This table 
lists the precise gluings of tetrahedron faces. For instance, the second cell in the bottom row 
indicates that face 2 of tetrahedron 3 is glued to tetrahedron 2, in such a way that vertices 1, 3, 4 
of tetrahedron 3 map to vertices 4, 2, 3 of tetrahedron 2 respectively. This same gluing can be 
seen from the other direction by examining the first cell in the middle row. 
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Face 1 
Vertices 234 


Face 2 
Vertices 134 


Face 3 
Vertices 124 


Face 4 
Vertices 123 


Tet. 1 
Tet. 2 
Tet. 3 


Tet. 1 


231 


Tet. 2 
Tet. 1 


134 
134 


Tet. 3 
Tet. 2 
Tet. 1 


124 
123 
124 


Tet. 1 
Tet. 2 
Tet. 3 


423 
124 
234 


Tet. 3 


341 


Tet. 3 


123 


Tet. 2 


423 



Table 2: The tetrahedron face gluings for an example 3-tetrahedron triangulation 

It is simple to see that the labelling for this triangulation is canonical. The sequence 
^1,1, . . . , An^4 is 1,2, 3, 1, 3, 1, 2, 2, 3, 2, 1, 3 (reading tetrahedron numbers from left to right and 
then top to bottom in the table), and tetrahedron 2 first appears before tetrahedron 3 as re- 
quired. Looking closer, the first appearance of tetrahedron 2 is in the second cell of the top row 
where vertices 1, 3, 4 map to 1, 3, 4, and the first appearance of tetrahedron 3 is in the following 
cell where vertices 1, 2, 4 map to 1, 2, 4. In both cases the gluings use the identity map. 

Lemma 5. For any triangulation T of size n, there are precisely 2An canonical lahellings ofT, 
and these can be enumerated in 0{n^ ^ogn) time. 

Proof. In summary, we can choose any of the n tetrahedra to label as tetrahedron 1, and we 
can choose any of the 4! = 24 labellings of its four vertices. From here the remaining labels are 
forced, and can be deduced in 0{n\ogn) time. The full proof is given in the appendix. □ 

Definition (Isomorphism signature). For any triangulation T of size n, enumerate all 24n 
canonical labellings of T, and for each canonical labelling encode the full set of face gluings as 
a sequence of bits. We define the isomorphism signature to be the lexicographically smallest of 
these 24n bit sequences, and we denote this by cr(T). 

To encode the full set of face gluings for a canonical labelling, we could simply convert a 
table of gluing data (such as Table [5]) into a sequence of bits. For practical implementations we 
use a more compact representation, which will be described in the full version of this paper. 

Theorem 6. Given two 3-manifold triangulations T and T' , we have a{T) — ct(7~') if and only 
if T and T' are isomorphic. 

Proof. It is clear that cr(T) — cr{T') implies that T and T' are isomorphic, since both signatures 
encode the same gluing data. Conversely, if T and T' are isomorphic then their 24n canonical 
labellings are the same (though they might be enumerated in a different order). In particular, 
the lexicographically smallest canonical labellings will be identical; that is, cr{T) = cr{T'). □ 

Theorem 7. Given a 3-manifold triangulation T of size n, the isomorphism signature cr[T) 
has 0{n\ogn) size and can be generated in 0{n^\ogn) time. 

Proof. To encode a full set of face gluings, at worst we require a table of gluing data such as 
Table [2l with An cells each containing four integers. Because some of these integers require 
O(logn) bits (the tetrahedron labels), it follows that the total size of ct{T) is 0(n log n). 

The algorithm to generate cr(T) is spelled out explicitly in its definition. The 24n canonical 
labellings of T can be enumerated in O(n^logn) time (Lemma [5]). Because a full set of face 
gluings has size 0{n\ogn), we can encode the 24n bit sequences and select the lexicographically 
smallest in 0(n^ ^ogn) time, giving a time complexity of Oin? \ogn) overall. □ 
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This space complexity of 0{n\ogn) is the best we can hope for, since Theorem 0] shows that 
the number of distinct isomorphism signatures for size n triangulations grows hke exp(0(n logn)). 

It follows from Theorems [5] and [7] that isomorphism signatures are ideal tools for constructing 
arcs in the Pachner graph, as explained at the beginning of this section. Moreover, the relevant 
definitions and results are easily extended to bounded and ideal triangulations (which are beyond 
the scope of this paper). We finish with a simple but important consequence of our results: 

Corollary 8. Given two 3-manifold triangulations T and T' each of size n, we can test whether 
T and T' are isomorphic in 0{n^ logn) time. 

4 Analysing the Pachner graph 

As discussed in the introduction, our focus is on one-vertex triangulations of the 3-sphere. We 
therefore direct our attention to ^i{S^), the restricted Pachner graph of the 3-sphere. 

In this section we develop algorithms to bound the shortest length and smallest excess height 
of any simplification path from a given node at level n of ^3^1(5*^). By running these algorithms 
over the full census of 31017 533 one- vertex triangulations of the 3-sphere (as described in 
Section [3. ip . we obtain a computer proof of the following results: 

Theorem 9. From any node at level n of the graph ^i{S'^) where 3 < n < 9, there is a 
simplification path of length < 13, and there is a simplification path of excess height < 2. 

The bound 3 < n is required because there are no simplification paths in ^i{S^) starting at 
level 2 or below (see Figure For n > 9 a computer proof becomes computationally infeasible. 

The results of Theorem[5]are astonishing, especially in light of Mijatovic's super-exponential 
bounds. Furthermore, whilst it can be shown that the excess height bound of < 2 is tight, 
the length estimate of < 13 is extremely rough: the precise figures could be much smaller still. 
These results have important implications, which we discuss later in Section O 

In this section we describe the algorithms behind Theorem[51 and we present the experimental 
results in more detail. Our algorithms are constrained by the following factors: 

• Their time and space complexities must be close to linear in the number of nodes that 
they examine, due to the sheer size of the census. 

• They cannot loop through all nodes in ^i{S^), since the graph is infinite. They cannot 
even loop through all nodes at level n > 10, since there are too many to enumerate. 

• They cannot follow arbitrary breadth-first or depth-first searches through ^i{S^), since 
the graph is infinite and can branch heavily in the upward direction0 

Because of these limiting factors, we cannot run through the census and directly measure 
the shortest length or smallest excess height of any simplification path from each node. Instead 
we develop fast, localised algorithms that allow us to bound these quantities from above. To our 
delight, these bounds turn out to be extremely effective in practice. The details are as follows. 

general, a node at level n can have up to 2n distinct neighbours at level (n + 1). 



9 



4.1 Bounding excess heights 



In this section we compute bounds _ff„ so that, from every node at level n of the graph ^i{S^), 
there is some simphfication path of excess height < Hn- As in Theorem |9l we compute these 
bounds for each n in the range 3 < n < 9. 

Algorithm 10 (Algorithm for computing This algorithm runs by progressively building a 

subgraph G C ^i{S^). At all times we keep track of the number of distinct components of G 
(which we denote by c) and the maximum level of any node in G (which we denote by £). 

1. Initialise G to all of level n of 3^i{S^). This means that G has no arcs, the number of 
components c is just the number of nodes at level n, and the maximum level is £ — n. 

2. While c > 1, expand the graph as follows: 

(a) Construct all arcs from nodes in G at level i to (possibly new) nodes in ^i(S^) at 
level £ + 1. Insert these arcs and their endpoints into G. 

(b) Update the number of components c, and increment £ by one. 

3. Once we have c — 1, output the final bound iJ„ — £ ~ n and terminate. 

In step 2(a) we construct arcs by performing 2-3 moves. We only construct arcs from nodes 
already in G, which means we only work with a small portion of level £ for each £ > n. In 
step 2(b) we use union- find to update the number of components in small time complexity. 

It is clear that Algorithm [TOl is correct for any n > 3: once we have c = 1 the subgraph G 
is connected, which means there is a path from any node at level n to any other node at level 
n. By Theorem [T] at least one such node allows a 3-2 move, and so any node at level n has a 
simplification path of excess height < £. 

However, it is not clear that Algorithm 1101 terminates: it might be that every simplification 
path from some node at level n passes through nodes that we never construct at higher levels 
£ > n. Happily it does terminate for all 3 < n < 9, giving an output of H2 = 2 each time. 
Table |3] shows how the number of components c changes throughout the algorithm in each case. 



Input level n 


3 


4 


5 


6 


7 


8 


9 


Value of c when £ = n 
Value of c when £ = n + 1 
Value of c when £ = n + 2 


20 
8 
1 


128 
50 
1 


1297 
196 
1 


13 660 
1074 
1 


169 077 
7 784 
1 


2 142 197 
64 528 
1 


28 691 150 
557 428 
1 


Final bound Hn 


2 


2 


2 


2 


2 


2 


2 



Table 3: Results obtained when running Algorithm [10] for 3 < n < 9 



It is straightforward to show that the space and time complexities of Algorithm [TU] are linear 
and log-linear respectively in the number of nodes in G (other small polynomial factors in n and 
£ also appear). Nevertheless, the memory requirements for n = 8 were found to be extremely 
large in practice (^29 GB), and for n = 9 they were too large for the algorithm to run (estimated 
at 400-500 GB). In the case of n = 9 a two-phase approach was necessary: 

1. Use Algorithm [TUl for the transition from level n to level n+ 1, and terminate if Hn ~ 1. 
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2. From each node v at level n + I, try all possible combinations of a 2-3 move followed by 
a 3-2 move. Let w be the endpoint of such a combination (so w is also a node at level 
n + I). If w € G then merge the components and decrement c if necessary. Otherwise do 
nothing (since w would never have been constructed in the original algorithm). 

3. If c = 1 after this procedure then output _ff„ — 2; otherwise terminate with no result. 

It is important to note that, if this two-phase approach does output a result, it will always 
be the same result as Algorithm [101 Essentially Step 2 simulates the transition from level n + 1 
to n -|- 2 in the original algorithm, with the advantage of a much smaller memory footprint (since 
it does not store any nodes at level n + 2), but with the disadvantage that it cannot move on to 
level n + 3 if required (and so it cannot output any result if iJ„ > 2). 

Of course by the time we reach n = 9 there are reasons to suspect that Hn ~ 2 (following 
the pattern for 3 < n < 8), and so this two-phase method seems a reasonable (and ultimately 
successful) approach. For n = 9 the memory consumption was '-^SOGB, which was (just) within 
the capabilities of the host machine. 

4.2 Bounding path lengths 

Our next task is to compute bounds L„ so that, from every node at level n of ^i{S^), there is 
some simplification path of length < Once again we compute L„ for 3 < n < 9. 

Because it is infeasible to perform arbitrary breadth-first searches through ^i{S^), we only 
consider paths that can be expressed as a series of jumps, where each jump involves a pair of 2-3 
moves followed by a pair of 3-2 moves. This keeps the search space and memory usage small: 
we always stay within levels n, n + 1 and n + 2, and we never need to explicitly store any nodes 
above level n. On the other hand, it means that our bounds L„ are very rough — there could be 
much shorter simplification paths that we do not detect. 

Algorithm 11 (Algorithm for computing -L„). First identify the set I of all nodes at level n 
of g^i{S^) that have an arc running down to level n — 1. Then conduct a breadth-first search 
across level n, beginning with the nodes in I and using jumps as the steps in this breadth-first 
search. If j is the maximum number of jumps required to reach any node in level n from the 
initial set I, then output the final bound Ln = 4j + 1. 

To identify the initial set I we simply attempt to perform 3-2 moves. When we process each 
node V, we must enumerate all jumps out from v; that is, all combinations of two 2-3 moves 
followed by two 3-2 moves. The number of such combinations is O(n^) in general. 

This time we can guarantee both correctness and termination if 3 < n < 9. Because n > 3 
the initial set / is non-empty (Theorem [1]) , and from our height experiments in Section l4.ll we 
know that our search will eventually reach all of level n. It follows that every node at level n of 
^i(S'^) has a path of length < 4j to some v G I, and therefore a simplification path of length 
< 4j -|- 1. Table 2] shows how the search progresses for each n. 

This time the space and time complexities are linear and log-linear respectively in the number 
of nodes at level n (again with further polynomial factors in n). This is considerably smaller 
than the number of nodes processed in Algorithm 1101 ^-i^d so for Algorithm 1111 memory is not a 
problem: the case n = 9 runs in under 4 GB. 
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Input level n 


3 


4 


5 


6 


7 


8 


9 


Size of I 




3 


46 


504 


6 975 


91283 


1 300 709 


18 361866 


Nodes remaining 




17 


82 


793 


6 685 


77 794 


841 488 


10 329 284 


Nodes remaining 


after 1 jump 


3 


1 


19 


75 


496 


4 222 


31250 


Nodes remaining 


after 2 jumps 








1 


1 





6 


12 


Nodes remaining 


after 3 jumps 























Final bound L„ 


9 


9 


13 


13 


9 


13 


13 



Table 4: Results obtained when running Algorithm [11] for 3 < n < 9 



4.3 Parallelisation and performance 

For n = 9, both Algorithms [10] and [TT] have lengthy running times: Algorithm [TOj requires a 
very large number of nodes to be processed at levels 9, 10 and 11 of the Pachner graph, and 
Algorithm 1111 spends significant time enumerating the 0(n*) available jumps from each node. 

We can parallelise both algorithms by processing nodes simultaneously (in step 2 of Algo- 
rithm[Tni and during each stage of the breadth-first search in Algorithm [TT|) . We must be careful 
however to serialise any updates to the graph. 

The experiments described here used an 8-core 2.93 GHz Intel Xeon X5570 CPU with 72 GB 
of RAM (using all cores in parallel). With the serialisation bottlenecks, Algorithms [TOl and 
[TT] achieved roughly 90.5% and 98.5% CPU utilisation for the largest case n = 9, and ran for 
approximately 6 and 15 days respectively. All code was written using the topological software 
package Regina [1,3. 

5 Discussion 

As we have already noted, the bounds obtained in Section[3]are astonishingly small. Although we 
only consider ri < 9, this is not a small sample: the census includes ^ 150 million triangulations 
including ^ 31 million one-vertex 3-spheres; moreover, nine tetrahedra are enough to build 
complex and interesting topological structures [U [13]. Our results lead us to the following 
conjectures: 

Conjecture 1. From any node at any level n > 3 of the graph 3^i(S'^) there is a simplification 
path of excess height < 2. 

If true, this result (combined with Theorem]!} would reduce Mijatovic's bound in TheoremJJ] 
from exp(0(ri^)) to exp(0(n log n)) for one- vertex triangulations of the 3-sphere. Furthermore, 
it would help explain why 3-sphere triangulations are so easy to simplify in practice. 

There are reasons to believe that a proof might be possible. As a starting point, a simple 
Euler characteristic argument shows that every closed 3-manifold triangulation has an edge of 
degree < 5; using at most two "nearby" 2-3 moves, this edge can be made degree three (the 
setting for a possible 3-2 simplification). The details will appear in the full version of this paper. 

Conjecture 2. From any node at any level n > 3 of the graph 3^i(S'^) there is a simplification 
path of length < 13. 

This is a bolder conjecture, since the length experiments are less consistent in their results. 
However, the fact remains that every 3-sphere triangulation of size n < 9 can be simplified after 
just three jumps, and this number does not rise between n — b and n = 9. 
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If true, this second conjecture would yield an immediate polynomial-time 3-spliere recogni- 
tion algorithm: for any triangulation of size n > 3 we can enumerate all 0{n^^^) combinations 
of three jumps, and test each resulting triangulation for a 3-2 move down to n — 1 tetrahe- 
dra. By repeating this process n — 2 times, we will achieve either a recognisable 2-tetrahedron 
triangulation of the 3-sphere, or else a proof that our input is not a 3-sphere triangulation. 

Even if Conjecture [His false and the length bounds do grow with n, this growth rate appears 
to be extremely slow. A growth rate of L„ S O(logn) or even 0(-\/n) would still yield the first 
known sub-exponential 3-sphere recognition algorithm (using the same procedure as above), 
which would be a significant theoretical breakthrough in algorithmic 3-manifold topology. 

Looking forward, it is natural to ask whether this behaviour extends beyond the 3-sphere 
to triangulations of arbitrary 3-manifolds. Initial experiments suggest "partially" : the Pachner 
graphs of other 3-manifolds also appear to be remarkably well-connected, though not enough to 
support results as strong as Conjectures [T] and [5] above. We explore these issues further in the 
full version of this paper. 
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Appendix: Additional proofs 

Here we offer full proofs for Theorem 2] and Lemma [5l which were omitted from the main text 
to simplify the exposition. 

Theorem 4. The number of distinct isomorphism classes of 3-manifold triangulations of size 
n grows at an asymptotic rate o/exp(0(nlogn)). 

Proof. An upper bound of exp(0(nlogn)) is easy to obtain. If we count all possible gluings of 
tetrahedron faces, without regard for isomorphism classes or other constraints (such as the need 
for the triangulation to represent a closed 3-manifold) , we obtain an upper bound of 

[(4n - 1) X (4n - 3) X • • • X 3 X 1] • 6^" < [Anf ■ 6^" G cxp(0(n log n)). 

Proving a lower bound of exp(J7(nlogn)) is more difficult — the main complication is that 
most pairwise identifications of tetrahedron faces do not yield a 3-manifold at all [5]. We work 
around this by first counting 2-manifold triangulations (which arc much easier to obtain), and 
then giving a construction that "fattens" these into 3-manifold triangulations without introduc- 
ing any unwanted isomorphisms. 
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To create a 2-nianifold triangulation of size 2m (the size must always be even), we identify 
the 6to edges of 2m distinct triangles in pairs. Any such identification will always yield a closed 
2-manifold (that is, nothing can "go wrong", in contrast to the three-dimensional case). 

There is, however, the issue of connectedness to deal with (recall from the beginning of 
Section [5] that all triangulations in this paper are assumed to be connected). To ensure that a 
labelled 2-manifold triangulation is connected, we insist that for each fc = 2, 3, . . . , 2m, the first 
edge of the triangle labelled k is identified with some edge from one of the triangles labelled 
l,2,...,fc — 1. Of course many connected labelled 2-manifold triangulations do not have this 
property, but since we are proving a lower bound this does not matter. 

We can now place a lower bound on the number of labelled 2-manifold triangulations. First 
we choose which edges to pair with the first edges from triangles 2,3,..., 2m; from the property 
above we have 3 x 4 x . . . x 2m x (2m -I- 1) = ^ (2m -f 1) ! choices. We then pair off the remaining 
2m -I- 2 edges, with (2m + 1) x (2m — l)x...x3xl = (2m -I- l)!/2™m! possibilities overall. 
Finally we note that each of the 3m pairs of edges can be identified using one of two possible 
orientations. The total number of labelled 2-manifold triangulations is therefore at least 

(2m -hi)! (2m + 1)! _ (2m + 1)! • (2m + 1)! • 22™ 

2 2™m! ' " 2 • ml ' 

Each isomorphism class can contain at most (2m)! • 6^™ labelled triangulations, and so the 
number of distinct isomorphism classes of 2-manifold triangulations is bounded below by 

(2m + 1)! • (2m + 1)! • 2^™ _ (2m + 1) • (2m + 1)! 
2 • to! • (2m)! • 6^™ ^ 2 • m! • 3^" 

> (2m + 1) X 2m X • • • X (to + 2) X (to + 1) X (5)™ 

>(TO+l)"+l.(i)" 

e exp(51(TOlogTO)). 

We fatten each 2-manifold triangulation into a 3-manifold triangulation as follows. Let F 
denote the closed 2-manifold described by the original triangulation. 

1. Replace each triangle with a prism and glue the vertical faces of adjacent prisms together. 



as illustrated in Figure 4(a) This represents a bounded 3-manifold, which is the product 
space F X I. 

2. Cap each prism at both ends with a triangular pillow, as illustrated in Figure [4(b)[ The 
two faces of each pillow are glued to the top and bottom of the corresponding prism, 
effectively converting each prism into a solid torus. This produces the closed 3-manifold 



F X S-^, and the complete construction is illustrated in Figure 4(c) 



3. Triangulate each pillow using two tetrahedra, which are joined along three internal faces 
surrounding an internal vertex. Triangulate each prism using 14 tetrahedra, which again 
all meet at an internal vertex. Both triangulations are illustrated in Figure [4(d) [ 

If the original 2-manifold triangulation uses 2to triangles, the resulting 3-manifold triangu- 
lation uses n = 32to tetrahedra. Moreover, if two 3-manifold triangulations obtained using this 
construction are isomorphic, the original 2-manifold triangulations must also be isomorphic. 
The reason for this is as follows: 
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(c) The complete construction (d) Triangulating prisms and pillows 

Figure 4: Fattening a 2-nianifold triangulation into a 3-nianifold triangulation 



• Any isomorphism between two such 3-manifold triangulations must map triangular pillows 

to triangular pillows. This is because the internal vertex of each triangular pillow meets 
only two tctrahcdra, and no other vertices under our construction have this property. 

• By "flattening" the triangular pillows into 2-dimensional triangles, we thereby obtain an 
isomorphism between the underlying 2-manifold triangulations. 

It follows that, for n = 32m, we obtain a family of exp(il(mlogm)) = exp{il{nlogn)) 
pairwise non-isomorphic 3-manifold triangulations. 

This result is easily extended to n ^ mod 32. Let Vn denote the number of distinct 
isomorphism classes of 3-manifold triangulations of size n. 

• Each triangulation of size n has at least n — 1 distinct 2-3 moves available (since any face 
joining two distinct tetrahedra defines a 2-3 move, and there are at least n—1 such faces). 

• On the other hand, each triangulation of size n+1 has at most 6(n-|- 1) distinct 3-2 moves 

available (since each 3-2 move is defined by an edge that meets three distinct tetrahedra, 
and the triangulation has at most 6(n + 1) edges in total). 

It follows that Vn+i > Vn ■ g^^^ > Vn/18 for any n > 1. This gives ^32™+*: > ^32^/18^^ 
for sufficiently large m and all < A; < 32, and so we obtain Vn € exp(rj(nlogn)) with no 
restrictions on n. □ 

Remark. Of course, we expect that Vi+i » Vn (and indeed we see this in the census). The 
bounds that we use to show Vn+i > in the proof above are very loose, but they are 

sufficient for the asymptotic result that we seek. 

Lemma 5. For any triangulation T of size n, there are precisely 24n canonical labellings ofT, 
and these can he enumerated in 0(n^ log n) time. 
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Proof. For n = 1 the result is trivial, since all 24 = 4! possible labellings are canonical. For 

n > 1 wc observe that, if we choose (i) any one of the n tetrahedra to label as tetrahedron 1, 
and (ii) any one of the 24 possible labellings of its four vertices, then there is one and only one 
way to extend these choices to a canonical labelling of T. 

To see this, we can walk through the list of faces -Fi.ij -^1,2, ^i,3) Fi,4, F2,i, ■ ■ ■ , Fn,4, where 
Ft^i represents face i of tetrahedron t. The first face amongst -Fi^i, . . . , Fi^^ that is joined to 
an unlabelled tetrahedron must in fact be joined to tetrahedron 2 using the identity map; this 
allows us to deduce tetrahedron 2 as well as the labels of its four vertices. 

We inductively extend the labelling in this manner: once we have labelled tetrahedra 1,. . . ,k 
and their corresponding vertices, the first face amongst Fi^i, . . . , Fk^4 that is joined to an unla- 
belled tetrahedron must give us tetrahedron k + 1 and the labels for its four vertices (again using 
the identity map). The resulting labelling is canonical, and all of the labels can be deduced in 
0(n log n) time using a single pass through the list Fi^i, . . . , F„^4. The logn factor is required 
for manipulating tetrahedron labels, each of which requires O(logn) bits. 

It follows that there are precisely 24n canonical labellings of T, and that these can be 
enumerated in 0{n^ logn) time using 24n iterations of the procedure described above. □ 
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